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The title compound is investigated by specific heat measurements in the normal and supercon- 
ducting state supplemented by upper critical field transport, susceptibility and magnetization mea- 
surements. Prom a detailed analysis including also full potential electronic structure calculations 
for the Permi surface sheets, Permi velocities and partial densities of states the presence of both 
strong electron-phonon interactions and considerable pair-breaking has been revealed. The spe- 
cific heat and the upper critical field data can be described to first approximation by an effective 
single band model close to the clean limit derived from a strongly coupled predominant hole sub- 
system with small Permi velocities. However, in order to account also for Hall-conductivity and 
thermopower data in the literature, an effective general two-band model is proposed. This two-band 
model provides a flexible enough frame to describe consistently all available data within a scenario of 
phonon mediated s-wave superconductivity somewhat suppressed by sizeable electron-paramagnon 
or electron-electron Coulomb interaction. Por quantitative details the relevance of soft phonons and 
of a van Hove type singularity in the electronic density of states near the Permi energy is suggested. 
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I. INTRODUCTION 

The recent discovery of superconductivity in the in- 
termetaUic antiperovskite compound MgCNi^i with a 
superconducting transition temperature of ~ 8 K is 
rather surprising considering its high Ni content. There- 
fore it is expected that this compound is near a ferromag- 
netic instabihty which might be reached by hole doping 
on the Mg sites The possibihty of unconventional super- 
conductivity due to the proximity of these two types of 
collective order has attracted great interest in the elec- 
tronic structure and the physics of the pairing mecha- 
nism. 

Band structure calculations^iMi^ for MgCNis re- 
vealed a domination of the electronic states at the 
Fermi surface by the 3d orbitals of Ni, suggest- 
ing presence of ferromagnetic spin fluctuations iSii 
^^C NMR measurements^ result in Fermi liquid be- 
havior with an electronic crossover at T w 50 K, 
growing formation of spin fluctuations below T w 
20 K. Resistivity measurementsr*2i£ measurements of 
the thermopower, the thermal conductivity and the 
magnetoresistancejiS doping experimentsSill and magne- 
tization measurementsii are consistent with this inter- 
pretation. 

MgCNig can be considered as the 3-dimensional ana- 
logue of the quaternary layered transition metal borocar- 
bides which exhibit superconducting transition tempera- 
tures up to Tc — 23 K. In spite of the much lower of 
MgCNig, its upper critical field Hc2 at low temperatures. 



i/c2(0) = 8...15 T 



8.12.13.14.15 



is comparable with that 



of the borocarbides or even higher. However, a rather 
different shape, especially near T^, for the temperature 
dependence of Hc2{T) is observed for these compounds. 
The Hc2{T) dependence of MgCNig is similar to that 
of usual (standard) intermetaUic superconductors which 
are described reasonably well within the isotropic single- 
band approximation and exhibit a steep slope of i?c2 (T) 
at Tc- 

Through analysis of specific heat data, MgCNig 
was characterized in the framework of a conventional, 
phonon-mediated pairing both as a moderatei*i^ and as 
a strongiiifi coupling superconductor. Strong coupHng 
is also suggested by measurements of the thermopoweniS 
and the large energy gap determined from tunneling 
experimentsji^ The question of the pairing symmetry is 
controversially discussed in the literature. ^^C NMR 
experimentsfi specific heat measurements!^ and tun- 
neling spectrciii support s-wave pairing in MgCNig, 
whereas earher tunneling spectrai^ and penetration 
depth measurementsi^ have been interpreted in terms of 
an unconventional pairing state. Recent measurements 
of the critical current of MgCNig may be interpreted in 
the latter sense, too^ 

In the present investigation, specific heat data of 
MgCNig in the normal and superconducting state were 
analyzed in detail with the aid of a realistic phonon 
model an d strong coupling corrections as suggested by 
ICarbottd ^ The results are brought into accordance with 
the two-band character of MgCNig emerging from band 
structure calculations and a parallel analysis of the up- 
per critical field i?c2(0), in order to find out a consistent 
physical picture explaining at least qualitatively various 
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FIG. 1: Partial density of states of the two bands in MgCNig 
corresponding to the two Fermi surface sheets shown in Fig. 
|2| Dotted line: total density of states. 
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FIG. 2: Fermi surface sheets of MgCNig. Fermi velocities 
are measured in different colors (see scales below the figure) 
in units of 10^ cm/s, i.e. blue color stands for slow and red 
color for fast quasiparticles. Upper panel and lower left panel: 
hole sheets corresponding to "band 1" in Fig. Q Lower right 
panel: electron sheets corresponding to "band 2" in Fig. U 
The right panels present the Fermi velocity distribution of 
the two sheets on the same absolute scale to demonstrate the 
slow (heavy) character of the holes. Yellow color: sides of 
filled electrons. The left panel shows the tip-distribution in 
the hole sheets on a smaller scale in more detail. 



available experimental results. 



II. ESSENTIALS OF THE THEORETICAL 
ELECTRONIC STRUCTURE 

Following previous work of one of the authors (H.R.)^ 
in the present section we remind the reader on some es- 
sential features and point out important, but neverthe- 



less so far unpublished details of the electronic structure 
of MgCNig which are crucial for a proper interpretation 
of the specific heat (total and Fermi surface sheet (FSS) 
related partial densities of states (DOS)), upper critical 
field and transport data (topology of the Fermi surface 
and the magnitude of the Fermi velocities) . Among var- 
ious band structure calculations there is general consen- 
sus about the qualitative topology of the Fermi surface 
and the presence of a strong peak (Van Hove Singular- 
ity) in N{E) slightly below the Fermi energy. At the 
same time there are clear differences with respect to the 
magnitudes of 7V(0) = 4.8 states/eV =11 mJ/molK^ 
(to be compared with 4.63 states/eV;^ 4.99 states/eV;^ 
5.34 states/eVifi) and especially with respect to the Stoner 
factor S = 3.3 (compared with 1.75^ 2.78^ to 5^) as 
well as to the distance of the DOS peak 42 meV (com- 
pared with 40 meV— to 80 meV— ) below Ep. The peak 
may be of relevance for a proper quantitative description 
of the electronic specific heat, transport data, magnetic 
properties, and the superconductivity. Last but not least, 
there is also a sizeable disparity on the magnitude of the 
electron-phonon coupling constant Aph (ranging between 
0.8 . . . 2.0) mainly caused by the poor knowledge of the 
phonon energies and possible lattice anharmonicitiesiS^ 

Our results have been obtained by a band structure 
calculation code using the full-potential nonorthogonal 
local-orbital (FPLO) minimum-basis schemeiSi There are 
about 0.285 charges per unit cell with exactly equal num- 
bers of holes and electrons, i.e. rih = rig] which follows 
from the even number of electrons per unit cell. In other 
words, MgCNi3 is a so called compensated metal which 
must be described in terms of multi-band model by defi- 
nition (per se). Thus, it makes sense to start with a two- 
band model. The generalization to any higher multi-band 
scenario is straightforward. A standard single band sys- 
tem with an even number of electrons per unit cell would 
be an insulator and not a superconductor. Thus metal- 
icity is achieved owing to the two-band character which 
leads to electron and hole derived Fermi surface sheets 
(FSS). The total DOS iV(0) at the Fermi level can be 
decomposed into a roughly 85% and a 15% contribution 
stemming from two hole and two electron sheets of the 
Fermi surface, respectively (see Figs. n],l2l)- On the Fermi 
surface sheets shown in Fig. |2the sides of filled electron 
states are shown in yellow/gold color. The two types 
of hole sheets are formed by eight droplets (ovoids) ori- 
ented along the spatial diagonals of the cube, i.e. along 
the F-R lines and six FSS with a "four-leaved clover"- 
like shape centered at the X-points in the middle of the 
faces of the cube (see Fig. [SJ. The coordinates of the 
symmetry points read F = (0,0,0), R = (0.5,0.5,0.5), 
X = (0.5, 0, 0) and M = (0.5, 0.5, 0) (all given in units of 
27r/a, where a = 0.381 nm is the lattice constant). The 
FSS with electron character are given by the rounded 
cube centered at T and 12 thin jungle gims spanning 
from R to M. 

The band structure calculations provide us directly 
with several material parameters (total and partial den- 
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sities of states, Fermi velocities, etc.) important for the 
understanding of superconductivity and electronic trans- 
port properties. For instance the transport properties 
are described by quadratically averaged Fermi velocities: 
(^^^)pgg whereas the upper critical field is described by 
averages of the type (l/w^)pgg which yields a smaller 
effective velocity in general. Using the general defini- 
tions of the local density of states (in k -space) and 
those of m*'' and the first moments of the Fermi velocity 



respectively, we have 
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where i = el, h, tot and Sf^i denotes the area of the i*^ 
Fermi surface sheet and the effective quantity is related to 
the linearly averaged value v adopting a simple estimate 
as 



Vi J dSi ~v'i(^k^ 



Vhc2 



SF,i 

V 1 — {5v/v)^ 



^hc2 'j' 



where 5v is the halfwidth of the v-p distribution. For 
the two above mentioned subgroups of quasiparticles we 
estimate: V2 = 3.9 x lO'' cm/s and vi = 1.2 x 10'' cm/s, 
Dtr,i = 1.11 X 10'' cm/s, Uhc2,i = 1-07 X 10'^ cm/s and 
Vtr:,2 = 4.89 X 10^ cm/s, where 5vi = 4 x 10^ cm/s and 
6v2 = 1.1 X 10'' cm/s have been used (compare also Fig. 

El. 

Finally, in the isotropic single band (ISB) model real- 
ized in the extreme dirty limit of superconductivity one 
arrives at 



ISB 



which yields Wtr.iSB = 2.15 x 10*" cm/s in accordance 
with Ref. 0. The corresponding plasma energy amounts 
hio-p\ = 3.17 eV close to 3.25 eV given in Ref. 0. Natu- 
rally, the total plasma frequency Wpi can be also decom- 
posed into the plasma frequencies of both subsystems 



pi 



V,2 



-'pl.l 



Thus we estimate Swpi^i = Swpi^h ~ 1-89 . . . 1.94 eV and 
tujjp\_2 = fiwpi,ei ~ 2.55 . . . 2.61 eV. From these partial 
plasma energies a useful relation between the scattering 
rates 7i and the conductivities <t\ (with i = 1,2) in both 
subsystems can be obtained: 



1^ ^ 1.816^ ^ 1.816^, 
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where p\ denotes the corresponding resistivity. In the 
present case the disorder is expected to be caused mainly 
by Mg and C related defects such as vacancies and inter- 
stitials. Therefore the ratio of the scattering rates might 
scale with the ratio of the non-Ni derived Mg and C or- 
bital partial densities of states at the Fermi level and the 
corresponding Fermi velocities. As a result we estimate 
from our LDA-FPLO calculations 



V 7h / LDA 



^el,Mg,c(0)«F,el 
^h,Mg,c(0)wF,h 



! 4.81. 



(1) 



Within this approach the corresponding mean free paths 
differ by a factor of 0.917 and a conductivity ratio of 
fh/cei = 1.403 would be expected. 

In the following analysis we usually make use of ?i = 
= A'o = 1 for the sake of simplification. 



III. EXPERIMENTAL 

Polycrystalline samples of MgCNia have been prepared 
by solid state reaction. In order to obtain samples with 
high Tc, we used an excess of carbon as proposed in 
Ref. H To cover the high volatility of Mg during sin- 
tering of the samples an excess of Mg is usedi In this 
study, a sample with the nominal formula Mgj^ 2Ci.6Ni3 
has been investigated and is denoted as MgC^ gNis. To 
prepare the sample a mixture of Mg, C and Ni powders 
was pressed into a pellet. The pellet was wrapped in a 
Ta foil and sealed in a quartz ampoule containing an Ar 
atmosphere at 180 mbar. The sample was sintered for 
half an hour at 600° C followed by one hour at 900° C. 
After a cooling process the sample was reground. This 
procedure was repeated two times in order to lower a 
possible impurity phase content. The obtained sample 
was investigated by x-ray diffractometry to estimate its 
quality. The diffractometer pattern (Fig. ^ shows small 
impurity concentrations mainly resulting from MgO and 
unreacted carbon crystallized as graphite (« 10 voL- 
%). The lattice constant of the prepared sample was 
determined to be a = 0.38107(1) nm using the Rietveld 
code FULLPROFi^ This indicates that the nearly single- 
phase sample corresponds to the superconducting mod- 
ification of MgC^NiaiS^ The superconducting transition 
of the sample was investigated by measurements of elec- 
trical resistance, ac susceptibility and specific heat. For 
the electrical resistance measurement a piece cut from 
the initially prepared pellet with 5 mm in length and a 
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FIG. 3: Rietveld refinement for the MgCj gNis sample. The 
crosses correspond to the experimental data. The black line 
shows the calculated pattern. The vertical bars give the Bragg 
positions for the main phase MgCNi3, for graphite and MgO 
(from top to bottom). The black line at the bottom of the plot 
gives the difference between the experimental and calculated 
pattern. 
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FIG. 4: Resistivity as a function of temperature of the 
MgC]^ ^^13 sample up to room temperature. The inset shows 
the superconducting transition region. 



B. Specific heat 



cross section of approximately 1 mm^ was measured in 
magnetic fields up to 16 T using the standard four probe 
method with current densities between 0.2 and 1 A/cm^. 
The ac susceptibility and the specific heat measurements 
were performed on other pieces from the same pellet in 
magnetic fields up to 9 T. 



IV. RESULTS 

A. Resistivity 

The temperature dependence of the electrical resis- 
tance of the investigated sample is shown in Fig. 31 A su- 
perconducting transition with an onset (midpoint) value 
of Tc = 7.0 K (6.9 K) is observed (see inset of Fig. 
which coincides with the onset of the superconducting 
transition of Tc = 7.0 K determined from ac susceptibil- 
ity. Its residual resistance ratio p(300K)/p(8K) = 1.85 
and the shape of the p{T) curve are typical for MgCNia 
powder samplespiiSi^ It should be noted, that the sample 
of Fig. 01 has a resistivity of psooK =2.1 mficm which is 
much too large in order to be intrinsic. 

A natural explanation for the high resistivity of the 
investigated sample which was not subjected to high 
pressure sintering is a relatively large resistance of 
the grain boundariesi^ In this context the recent low- 
resistivity th in film data (with po down to 20 pflcm) by 
lYoung et alife i are of interest, since in that case similar 
values of the upper critical field and just as in the best 
powder samples^ have been observed. 



Specific heat measurements were performed in order 
to get information about the superconducting transition, 
the upper critical field and the superconducting pairing 
symmetry and the strength of the electron-phonon cou- 
pling from thermodynamic data. In Fig. specific heat 
data, Cp/T vs. T^, are shown for applied magnetic fields 
up to 8 T. The previously mentioned (see Sec. Illljl 10 
vol.-% graphite imp urity contribution was subtracted ac- 
cording to Ref. |2a. 

The specific heat can be considered as a sum of a lat- 
tice contribution and a linear-in-T term which gives the 
electronic contribution with 7^ as the Sommerfeld pa- 
rameter: 

Qattice 

(T). (2) 

To extract the lattice contribution of the normal state 
specific heat the low temperature limit 

Clattice(r) = /3T^ (3) 

of the Debye model is usually applied. A fit of Eq. Q to 
the data is shown in Fig. [^resulting the parameters f3 ~ 
0.39 mJK^/mol and 7^ = 27.0 mJ/molK^. Notice, that 
the Sommerfeld parameter is connected to the electron- 
phonon coupling strength by 7;!^ = 7o (1 + '^ph)- With 
70 = 11 mJ/molK^ (Sec. m, one obtains Aph = 1.45 
in contradiction with recently reported medium coupling 
resultsiiii^ From the lattice contribution the Debye tem- 
perature 6q = 292 K was derived. Both parameters 
(6q, j^) are consistent with what has been reported so 
farpi^ii^ The fit describes the normal state data above 
Tc but its extrapolation to T = K obviously underesti- 
mates the high-field data (see Fig. 

The transition temperature, = 6.8 K, calculated 
from entropy conservation criterion, agrees well with the 
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FIG. 5: Specific heat data Cp/T vs. of MgCi.eNis mea- 
sured at various magnetic fields up to 8 T. The dashed line 
is a fit of the Debye approximation to the data for H = T 
above Tc- Its intersection with the Cp/T-axis gives the Som- 
merfeld parameter 7jj = 27 mJ/molK^ (see text). 



transition temperatures ~ 6.9 K and T^, = 7.0 K de- 
rived from resistance and from ac susceptibility data, re- 
spectively. 

The jump Ac(r = T^) of the specific heat is 
given by the difference between the experimental data, 
Cp(T) and the normal state specific heat contribution 
Cn(T). Notice that the experimental value of the jump, 
Ac(Tc)/ (7n?c) = 2.09 (derived from an entropy conserv- 
ing construction - see Sec. IV is strongly enhanced 
compared to the BCS value (1.43), indicating strong 
electron-phonon coupling. 



C. Superconducting transition and upper critical 
field 

The field dependence of the electrical resistance of our 
investigated sample is shown in Fig. El for several tem- 
peratures between 1.9 . . . 6.0 K. A sharp transition is ob- 
served. It remains sharp down to low temperatures. In 
Fig. the field values Hio, H^o and Hgo defined at 10 %, 
50 % and 90 % of the normal state resistance are plotted 
as function of temperature. Identical results have been 
found from resistance-vs. -temperature transition curves 
measured at different magnetic fields. Additionally, Fig. 
shows upper critical field data determined from ac sus- 
ceptibility measurements, H^^^, determined by an onset 
criterion. 

It is clearly seen that for the investigated sample i?™** 
agrees approximately with Hiq. A similar behavior was 
already observed for MgB2, whereas in the case of rare- 
earth nickel borocarbides the onset of superconductivity 
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FIG. 6: Resistivity of the MgCj^ gNia sample as a function of 
the applied magnetic field for various fixed temperatures as 
labeled. 



determined from ac susceptibility was typically found to 
agree well with the midpoint value (i/so) of the normal 
state resistivity. The width AH — Hqq — Hiq of the su- 
perconducting transition curves in Fig. M (and Fig. ^ 
remains, with AH ~ 0.6 T, almost unchanged down to 
low temperatures. A non-textured polycrystalline sample 
of a strongly anisotropic superconductor shows a gradual 
broadening of the superconducting transition with de- 
creasing temperature as was observed, for example, for 
MgBji^ Therefore, the nearly constant transition width 
AH observed for the investigated sample can be con- 
sidered as an indication of a rather small anisotropy of 
H,2{T) in MgCNig. 

The upper critical fields, Hc2(r), determined from the 
specific heat data, are shown in Fig. |H1 The Hc2{T) 
data obtained from the specific heat are located in the 
small field range between the Hg(j{T) and Hiq(T) curves 
determined from resistivity measurements (see Fig. ^ . 

The extrapolation of HgoiT) to T = K yields an up- 
per critical field of i?c2(0) ^ 11.0 T (see Fig. 01. The 
observed temperature dependence of the upper critical 
field is typical for Hc2{T) data reported for MgCNig so 
faiii^ and was describedSiiii^ within the standard WHH 
model2£ by conventional superconductivity. However, 
a quantitative analysis of Hc2 data presented in Sec. 
IVI Al shows that the magnitude of the upper critical field 
Hc2{0) at T = K can be understood only if strong 
electron-phonon coupling is taken into account. 



V. ANALYSIS 

A. Resistivity in the normal state 

The non-intrinsic origin of the residual resistivity 
PoK ~ 1.13 mficm follows from physically reasonable val- 
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FIG. 7: The upper critical field as a function of temperature. 
The circles show the midpoint of the transition (Hso). The 
two lines labeled Hio and Hgo denote 10 % and 90 % of normal 
state resistivity. The triangles represent the upper critical 
field from susceptibility measurements (onset values). The 
dashed line illustrates the extrapolation of the resistivity data 
to T = K. 
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FIG. 8: Comparison of upper critical field data determined 
from specific heat (•) and resistance measurements. Hio and 
Hqo were determined at 10% and 90% of the normal state 
resistivity, respectively. An entropy conserving construction 
was used to determine the upper critical field from the specific 
heat data of Fig. |^ 



ues for the mean free path I given by 

47rwF 



imp 9 



or in more convenient practical units 

vp [lO"^ cm/s] 



^imp[nm] = 4.9 x 10^ 



(fiWpi[eV]) poifi^cm] 



(4) 



TABLE I: Debye and Einstein temperatures with correspond- 
ing occupation numbers. T)i denote the acoustic phonons and 
Ei the optical phonons. O gives the corresponding tempera- 
ture and Ui is the grouping parameter, giving the number of 
modes found to have the same temperature. 



acoustic modes 
Dl D2 



El 



optical modes 
E2 E3 E4 



E5 



[K] 



129 
1 



316 
2 



86 163 
0.33 2.67 



256 
3 



472 
3 



661 
3 



with po as the residual resistivity. In the isotropic sin- 
gle band (ISB) case, with vp = WtrjSB ~ 2 x 10'' cm/s 
and the plasma energy fwjp\ « 3 eV (see Sec. one ar- 
rives at the conclusion, that in fact, in the mficm range, 
typical for most powder samples considered so far in 
the literature, even at T = K the mean free path I 
- given by Eq. ^ - would be smaller than the lattice 
constant a « 0.38 nm in obvious conflict with the well- 
known Joffe-Regel limit I > a^^^ In other words the 
maximal intrinsic resistivity is given by p™^'' [mficm] « 
1.3uF [10^ cm/s] / (fiwpi [eV])^ = 0.29 mQcm in the ISB 
dirty limit. Such resistivities (0.33 milcm) have been re- 
ported in mo st heavily (neutron) irradiated samples by 
iKarkin et alt i 



B. Specific heat in the normal state 

In order to describe the speciflc heat data in the normal 
state in an extended temperature range Tc < T < 30 K, 
the Debye low temperature limit approximation for the 
lattice contribution (see Eq. (pj) was replaced by 



Clatticel?") = CdCT) + Ce(T). 



Here, 



Cd(T) = s3^i 3i? 



T 



3 i-Bu./T 



dx- 



e X 



(e- - ly 



stands for the Debye model22i2i describing the 3 acoustic 
phonon branches, whereas the Einstein modelSSiSi 



CE(r) 



R 



9k 
T 



cxp(eE»/r) 
[exp(eE./r)-i]' 



describes the 12 optical branches. 

We found that the 9 energetically lowest phonons (3 
acoustic and 6 optical modes) are sufficient to describe 
the normal state specific heat up to T = 30 K. The result 
of the fit is shown in Fig. |2l The Sommerfeld parameter 
converged to = 31.4 mJ/molK^, greater than deter- 
mined from Fig. [3 Specific heat measurements up to 
T = 300 K on another piece from the initially prepared 
sample (which are not presented here) give the remain- 
ing 6 optical mode temperatures. The obtained Debye 
and Einstein temperatures and the belonging grouping 
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FIG. 9: Specific heat data Cp{T)/T vs. for zero magnetic 
field in the temperature range up to 30 K. The solid line is 
a fit of the lattice model (see text for details), showing very 
good agreement with the data for Tc < T < 30 K. 
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FIG. 10: Schematic phonon model spectrum -Fph(<j-') for 
MgCi gNis derived from fit parameters according to Tab. Q] 
The peak width of the optical modes was chosen arbitrarily 
as af — 1. 



parameters Vi are summarized in Tab. The phonon en- 
ergies are in good agreement with recent calculations 

The corresponding phonon spectrum has the formSi 



^^ph(^) 



^Dl H 1^1)2 



"di 



"D2 



E 



: exp 



{UJ - ilEiY 

2a? 



where 0{x) is the well known step-function and fl denotes 
the corresponding cut-off temperatures in meV.The re- 
sult including higher optical modes is shown in Fig. ^1 
Our model parameters even reproduce the rather com- 
plex phonon dispersion along the F-X direction in the 
first Brillouin zone at low phonon energies, as can be seen 
from Fig. IllL where th e used model is compared with cal- 
culations reported by llgnatov et aliB ^ Even though our 
model only involves constant and linear dispersion by 
the Einstein and Debye model, respectively, the calcu- 
lated phonon dispersion (right panel) is well reproduced 
(left panel), by means of superpositions of acoustic and 
optic phonon modes. The high-energy optic phonons ob- 
tained from the model are shifted to lower energies than 
predicted by the calculations. The shift is most proba- 
bly caused by anharmonic effects, which usually increase 
specific heat data at higher temperatures 

To investigate the electron-phonon coupHng strength, 
the electron-p honon interaction function a^Fph(a;) is 
of interest. The coupling function a^(w) is usually ex- 
tracted from tunneling measurements. In the case of A15 
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FIG. 11: Phonon dispersion along F-X direction in the first 
Brillouin zone, a) Our model with acoustic phonons (black 
lines) and optic phonons (dashed lines) according to Tab. U 
b) Calculations reported bv.I gnatov et all after Ref. [S^- 



Within this approach the logarithmically averaged mean 
phonon frequency win was determined from 



Win 



exp [ - — / duj — — — —Inu) 



2 / duj 

'0 



a^{u))F{uj) 



(5) 



compoundaSi and some borocarbidesjSi a^iu) is found to 
be of the form a^ioj) = 5/ ^/oj, with a scaling parameter 5. as win = 143 K. This frequency is used in the well known 
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McMillan formula (refined bv lAllen and Dvnesh ^ 



1.2 ^ 



1 + A 



ph 



Aph -M*(l + 0.6Aph)J 



(6) 



to estimate the electron-phonon coupling constant Aph- 
fi* denotes the usually weak Coulomb pseudopotential 
which has been chosen as ^i* — 0.13 in this case. 
With Tc = 6.8 K the electron-phonon coupling constant 
amounts Aph = 0.84, suggesting moderate coupling as 
proposed, for instance in Refs. Q and However, the 
low value of Aph estimated from Eq. ^ is in contradic- 
tion with our specific heat data as already mentioned in 
Sec. IIV Bl In particular, Aph = 1-45 was derived from 
the ratio 7j^/7o and also the high value of the super- 
conducting jump Ac(T(,)/(7Nr(,) = 2.09 indicates strong 
electron-phonon coupling. Strong electro n-phonon cou- 
pling was also predicted bv llgnatov et aliE ^ (Aph — 1.51) 
on the base of the calculations mentioned above. 

In this context a more precise analysis of the low tem- 
perate normal state specific heat data is required. As can 
be seen from the dashed line in Fig. El the extended 
lattice model does not describe the magnetic field data. 
Even larger deviations are observed if the experimental 
data are described within the lo w tempera ture limit of 
the Debye model (see Fig. iLin et al£ ^ who found 
a similar upturn of the experimental data at low tem- 
peratures tried to explain this behavior by the presence 
of Ni impurities. However, our x-ray analysis (see Fig. 
EJ shows no indicatio n for Ni impurities in our sample. 
Recently, IShan et alifc ^ found that the mentioned upturn 
can be easily reduced by lowering the carbon content. 
They attributed the observed upturn to some kind of bo- 
son mediated electron-electron interactions in MgC^Nis. 
This argument motivated us to search for other possi- 
ble sources to explain the low temperature upturn of the 
normal state specific heat data. 

The easiest explanation is an additional electron-boson 
interaction which may be an 

1. electron-phonon interaction originating from addi- 
tional phonon- softening of th e lowest acoustic mode 
(suggested hv llgnatov et aliE ^ and verified experi- 
mentally bv lHeid et aliE ^l and / or 

2. electron-paramagnon interaction (see Sec. HJ. 

Specific heat measurements let not clearly distinguish 
between these possible origins, but since magnetization 
measurements on our sample (not presented here) show 
increasing spin fiuctuations below ~ 30 K in accord with 
previous statements (see Sec. DJ, the focus in this pa- 
per lies on the electron-paramagnon interaction scenario. 
This is additionally supported by a small magnetic field 
dependence of the specific heat data, typically found in 
the presence of ferromagnetic spin fiuctuations. 

Within Eliashberg theory the renormalized normal 
state specific heat is described by the temperature de- 
pendent thermal mass Am* (T) /rriband- Its contribution 
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FIG. 12: Low temperature normal state total and electronic 
specific heat including field measurements (0.5 T (t), 2 T (★), 
4 T (■), 6 T (a) and 8 T (•)). Left panel: Specific heat data 
Cp(r)/r vs. T^. Dotted line: Extended lattice model de- 
scribing the zero field data (see Fig. O. Solid line: Fit of the 
model (including lattice and paramagnon contribution) to the 
data. Right panel: electronic specific heat 7N(r) vs. T in the 
normal state. Black line: Sommerfeld parameter 7N(r) of 
the model (see text for details), describing the observed up- 
turn of the specific heat at low temperatures. Dotted line: 
Qualitative model for spin fluctuations according to Eq. ((HJ. 



to the specific heat is given by 



A frr^ Am* (T) 
A7sf(r) = — 70, 



"^band 



with^ 



Am* (T) 

Wband 



6 



duja^F{uj) {-z 



2z^Im [tp'iiz)] - z^Re Wiiz)]} , 



where ip{iz) is the digamma function and z = 
uj/ (27rfcBT). The additional electron-paramagnon inter- 
action function is of the form 



a^Fgf = auj9 [flp 



(7) 



A corresponding fit results in a paramagnon-model tem- 
perature of flp ~ 2.15 meV =4> 9p 25 K with a thermal 
mass of Am*(r = K)/m « 0.43, which is of t he same 
order of magnitude as determined bv lSha.n et alfe Since 
this low energy excitation concerns the electronic part of 
the specific heat, we add it to the Sommerfeld parameter 
which than becomes temperature dependent. The elec- 
tronic contribution to the specific heat increases from ini- 
tially 7i^ = 31.4 mJ/molK^ to 7n (0) = 36.0 mJ/molK^ 
This is understandable since the paramagnon interaction 
dominates in the temperature range below 10 K. The 
black line in Fig. El shows the good agreement of this 
extended model with the experimental data in the low 
temperature region. The magnetic field dependence of 



9 



the paramagnons (which in addition can be temperature- 
dependent) is not included in the model. In the following 
7n (T) is denoted as 7n for the sake of simplicity. 
The usually applied model 



7N oc ST^ In (T/To) , 



to describe spin fluctuation behavior is shown in the right 
panel of Fig. El for comparison. 

At this point the question of the strength of the 
coupling may be rechecked. Including the additional 
electron-paramagnon interaction, the Allen-Dynes for- 
mula Eq. becomes 



Win 

1.2 



exp 



1 + A 



Aph - Agf - ^* (1 + 0.6Aph) 



(9) 



with A = Aph + Agf. Using win = 143 K, — 6.8 K, 
/i* = 0.13 and Agf = 0.43, the electron-phonon cou- 
pling constant rises to Aph — 1.85. Using this value, 
the electron-phonon interaction function based on the ap- 
proach a{u!) — S/^/uj can now be determined by scaling 
the factor S according to Eq. ij^l. The electron-boson 
interaction functions a^J^ph(t^) and a'^Fsf{uj) are shown 
in Fig. US 

The reliability of the model approach for the electron- 
phonon coupling function a{uj) = 8l\/u} can directly be 
checked from the band structure, using the ratio between 
the Sommerfeld parameter 7n(0) — 36.0 mJ/molK^ and 



the free electron parameter 
11.0 mJ/molK^ 



70 = i:'k'^N{E^)li 



7n(0) 
70 



(1 + Aph + Agf) . 



(10) 



With Agf K. 0.43, the electron-phonon coupling constant 
becomes Aph ~ 1.9, showing very good agreement be- 
tween both approaches. 

In the next section, the analysis of the specific heat in 
the normal state will be extended to the superconducting 
state. 



C. Specific heat in the superconducting state 

Fig. El shows the superconducting par t of the elec- 
tronic specific heat Ac(T) = Cp(r) — Cn(T), obtained 
from the zero-field data. The superconducting transition 
temperature = 6.8 K has been estimated by an en- 
tropy conserving construction (dashed line in Fig. I14II . 
This value agrees well with the transition temperatures 
Tc = 6.9 K and — 7.0 K, derived from resistance and 
from ac susceptibility data, respectively. The conserva- 
tion of entropy, AS'(r) = /g^°(Ac/T)dT, is shown in the 
inset of Fig. It was already mentioned, that the high 
value of the jump Ac(rc)/(7NT'c) ~ 2.09 found for the 
investigated sample can be explained by strong electron- 
phonon coupling. Nevertheless, we will start to analyze 
Ac(r) for T < T^jl within the BCS theory, since the de- 
viation from the weak-coupling temperature dependence 
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FIG. 13: Electron-phonon interaction function cP'F^\^(y>) for 
MgC]^ gNia. Phonon energies are marked by "Ei", respectively 
"Di" (see Fig. Ilf)l and Tab. Inset: Electron-paramagnon 
interaction function c?Fsi(y)) according to Eq. |J3, the para- 
magnon energy is marked by "P". 
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FIG. 14: Electronic specific heat data Ac/T vs. T in the 
superconducting state (filled circles). The solid line in the 
temperature range < T < 3.4 K corresponds to Eq. lllll . 
Dotted line: Entropy conserving construction to get the ide- 
alized jump. Inset: Entropy conservation for the electronic 
specific heat in the temperature range < T < Tc. 



of the gap is mainly restricted to the vicinity of the jump. 
The temperature dependence of Ac(T) = Cp(r) — Cn(r) 
in the weak-coupling BCS theory (Tc ^ win) is given by 
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FIG. 15: Normalized electronic specific heat contribution vs. 
Tc/T. The black line is a fit of Eq. Ijl2ll to the experimental 
data. The inset shows a fit of the two-band approximation 
given by Eq. l(T3)l. 



an approximative formula 



Ac(T) = 8.57nTc exp -0.82 



iBCS 



(0) 



TnT, (11) 



valid in the temperature range of 2 < T^/T < 6 corre- 
sponding in this case to 1 K < T < 3.4 K. Eq. lfTT|l 
can be fitted to the data by using the phenomenologi- 
cal gap 2Aexp/fcB?c = 3.75, slightly exceeding the BCS 
weak couphng prediction 2ABcs(0)/fcBTc — 3.52. The 
fit, which is shown as black fine in Fig. describes the 
experimental data in the range of2K<r<3.5K quite 
well. 

To examine the temperature dependence of the elec- 
tronic specific he at 



Cel 



cp(r) - 

Ac(T) 



^lattice 

(T) 



at = in detail, Cei(r)/7NTc is plotted logarithmically 
vs. Tc/T (Fig. El- The corresponding formula to Eq. 
ltTT|l reads 



Cel(r) 



8.5 exp 



_0.82^ 



(12) 



if Abcs(O) is replaced by Aexp- The black line in Fig. ^1 
is a fit of Eq. lfT2ll to the experimental data, which show 
an exponential temperature dependence at low tempera- 
tures (Tc/T > 2). This is a strong indication for s -wave 
superconductivity in MgC^ gNia. A natural explanation 
for the discrepancy from the expected BCS gap value 
of 2ABcs(0)/fcB7c = 3.52 to the experimentally found 
2Aexp/fcB7c = 3.75 emerges from 



1. the two-band character of MgC^ gNis and 

2. the enhanced electron-phonon coupling. 

The electronic specific heat can be analyzed within a two- 
gap model, simply by extending Eq. lfT2|l with a second 
gap 

Ai 



Cel(T) 



0.85 exp -0.82 



fcBT 



0.15 exp -0.82 



(13) 



using 85 % contribution for the hole band and 15 % con- 
tribution for the electron band. The dotted line in Fig. 
1 151 shows this fit for two gaps with 2Ai/fcBT'c = 3.67 (fit- 
ted parameter) and 2A2/fcB2c = 4.50 (fixed parameter). 
This situation is nearly indistinguishable from the single 
band model (see inset of Fig. El- The two-gap model 
is of considerable interest, since in this case the large 
gap found in recent tunne ling measurements of Ref. |lj 
(2A/fcBTc = 4.6) and Ref. [Tj (ZA/ZcbTc = 4.3) is repro- 
duced. 

Even the lower gap of the two-gap model exceeds the 
BCS weak-couphng limit. This and the strongly en- 
hanced specific heat jump Ac(T(,) are clear indications 
of strong electron-phonon coupling in accordance with 
our normal state specific heat analysis. Thus it is now 
straightforward to investigate the electron-phonon cou- 
pling strength and thus, the characteristic phonon fre- 
quency cjin, introduced in Sec. rVBl from the super- 
conducting state characteristics. The Eliashberg theory 
provides the following approximate formulas, which in- 
cludes strong couphng corrections within an isotropic sin- 
gle band model and links x = win/Tc to experimental 
thermodynamic quantitiestSS 

2A(0) 



Ac(r) 



kBTc 

Ac(re) 
Ac(re) 



- 7nT 



dT 



Tr 



3.53Bo {x) , 
1.43Bi (x) , 
-3.77^2 {x) , 
O.I68B3 (x) , 
0.576^4 (x) . 



(14a) 
(14b) 
(14c) 
(14d) 
(14e) 



The corresponding logarithmic correction terms are given 

by 



Bo{x) 


= 1 + 


12.5a;" 


-2 In-, 
2' 


(15a) 


Bi{x) 


= 1 + 


b2,x-^ 


In-, 

3' 


(15b) 


B2{,X) 


= 1 + 
= 1 - 


\YIx- 
I2.2x 


2ln— , 
2.9 


(15c) 
(15d) 


Bi{x) 


= 1- 


UAx 




(15e) 
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FIG. 16: Temperature dependence of the thermodynamic crit- 
ical field Hc{T) (filled circles) derived from the electronic spe- 
cific heat in the superconducting state using Eq. H16|l . Solid 
line (0 < T < 3.4 K): Single band model according to Eq. 
11211 . Dotted line: Idealized jump construction (see Fig. 11511 . 
Inset: Derivative dHc/dT vs. T (filled circles) and idealized 
jump (dotted line). 



Now, using Eq. H14a|l . = 6.8 K and the gap value 
Aexp(2 K) — 1.10 meV , one arrives at uji„ — 149 K. 

Using the value of the idealized jump of the specific 
heat, Ac(Tc)/ {-fj^Tc) = 2.09 in Eq. (tllbl with = 
6.8 K, Win = 88 K is derived. 

Comparing the linear slope of the idealized specific 
heat in the superconducting state of —6.7, obtained from 
Fig. IHwith Eq. iflicll . one gets win = 109 K. 

In view of strong-coupling effects the ratio 
1nT^/H^{Q), implying again only thermodynamic 
quantities, is of interest. Known superconductors show 
values between 0.17 . . .0.12 ranging from weak to strong 
coupling, respectively (see for example page 1086 of 
Ref. l2(i) . The thermodynamic critical field H^{T) can 
be determined with the help of the Gibbs free energy 
AF = -SAT - MdB as 



i?c(T) = V-SttAF. 



(16) 



AF is to be extracted from the specific heat in the su- 
perconducting state, Ac(r) = -Td'^{AF)/dT'^ . The 
temperature dependence of Hc{T) is shown in Fig. [161 
With i7c(0) = 179.6 mT we found jj^T^/H^{0) = 0.155. 
From Eq. l[T4d|l we get uji„ = 110 K (with = 6.8 K). 

Next, from the derivative of the thermodynamic 
critical field at zero temperature, dHc/dT the ratio 
ifc(O)/ {dHc{T)/dT\rp Tc) can be estimated. The value 
at T = T(. (of the idealized jump construction), amounts 
dHc{T)/dT\rp = 50.236 (see dashed line in the in- 
set of Fig. I16|l . Using the experimental value of 



Fc(0)/ (dfl"c(7')/dr|y^ Tc) = 0.525 in Eq. dill, a value 
of Win = 102 K is extracted. 

It should be noted that Eqs. Ill4a,ll and Ill4e|l can be 
used to estimate the value of the gap A(0) from the 
thermodynamic critical field i?c(0), due to similar de- 
pendences on strong coupHng corrections>2£ 



T dH,XT) 
HciO) dT 



T=T^ 



A(0) 
kv.Tr 



Using dH^{T)/dT\j.^ = 50.236, we get 2A(0)/fcBrc ~ 
3.80, agreeing well with the single band result 
2Aexp/fcBrc = 3.75 of Eq. ((T^. 

In summary, win was estimated from five different 
thermodynamic relations, only involving experimental 
results. The mean value w^T = (111 ± 23) K is in 
good agreement with calculations of llgnatov et alB ^ An 
overview of the results is given in Fig. El Note, that a 
similar analysis was already successfully used to describe 
some borocarbide superconductors i^^*^ 

The mean value wi^, derived from the superconduct- 
ing state is somewhat smaller than the normal state re- 
sult. Win = 143 K. This may be attributed to an ad- 
ditional phonon softening contribution or the approxi- 
mative approach of the electron-phonon coupling func- 
tion q;^(w) (see Sec. IVB|I . Nevertheless by checking Eq. 
(El with = (111 ± 23) K, Tc = 6.8 K, ^i* = 0.13 
and Agf = 0.43 the electron-phonon coupling constant 
becomes Aph = 1.9... 2. 3, whereas Aph = 1.85 was de- 
rived from Win = 143 K for the same parameters. It 
should be noted here, that Eqs. itTHl were derived assum- 
ing a small value for the Coulomb pseudopotential 
which is oversimplified considering enhanced electron- 
paramagnon coupHng found in this analysis. A rough cor- 
rection would shift the characteristic phonon frequency 
Win to slightly higher values and a coupling constant of 
''^ph ~ 1.9 seems to be most likely. Note that a similar 
analysis was already successfully used to describe some 
borocarbide superconductors i^^*^ 

The analysis of the thermodynamic properties of 
MgCNig presented so far clearly points to strong electron- 
phonon coupling. However, the temperature depen- 
dence of the thermodynamic critical field Hc{T) shown in 
Fig. El strongly deviates from analogous data for well- 
known strong coupHng superconductors such as Hg or 
Pb. Heir) is usually analyzed in terms of the deviation 
function D{t) = H^It) = H^iO) (l - t^) with t = T/T^. 
The deviation function of the above mentioned strong 
coupHng superconductors is positive and goes through 
a maximum at 0.5. The deviation function of 

MgCNig is shown in Fig. El Instead of the expected 
positive sign, D{t^) of MgCNig becomes negative already 
above about 0.3rc. The shape of the deviation function 
of MgCNig closely resembles that one of Nb having an 
electron-phonon coupling strength of Aph ~ 1.0. We re- 
mind the reader that the weak coupling BCS model yields 
a negative maximum deviation of « 3.8 % (see Fig. El- 
Thus, at first glance, our result seems to be in contradic- 
tion with the strong electron-phonon coupling suggested 
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FIG. 18: Deviation function of the thermodynamic critical 
field of MgC^Nis (filled circles) as function of {T jT^f . The 
solid line for < T < 0.34 K corresponds to Eq. Ill2|l . the dot- 
ted line corresponds to the idealized jump construction (see 
Fig. I14II . For comparison, the weak coupling BCS result*^ is 
shown. 



FIG. 17: Several thermodynamic quantities in dependence 
on the characteristic phonon frequency o^in according to Eqs. 
I|14|l and l|15|l . Filled circles: thermodynamic quantities es- 
timated for MgCi gNia from experimental data. Strong dis- 
crepancies are found within the low temperature Debye limit 
(filled squares). Note that the weak-coupling limit is reached 
in the asymptotic extrapolation — > oo. 



above. It turns out that this contradiction can be re- 
solved, taking into account a spHtting of the electron- 
phonon interaction function in a high and a low (soft) 
energy part. This is illustrated in Fig. where a two 
phonon peak spectrum with equal coupling strengths of 
both peaks located at cji and loi has been analyzed in the 
strong coupling case of Aph ~ 2 under consideration. The 
theoretical curves calculated within the ISB are shown 
for different frequency ratios LOijuj^- For ujiji^^ ~ 8, the 
"standard" strong coupling behavior, namely a positive 
deviation function, is completely removed and the devi- 
ation function becomes negative. Considering the low 
energy modes El and Dl, found in the analysis of the 
specific heat in the normal state (see Figs. 1101 and IT^ . 
this situation is easily imaginable to be valid in the case 

of MgCNig. 

In the superconducting state a Hnear-in-T electronic 
specific heat contribution ^{H)T arises from the normal 
conducting cores of the fiux lines for applied magnetic 
fields H> H,i. 

This contribution can be expressed as ^{H)T — 
Cp(T, 7J) — Cp(T, 0))^ where Cp(T, 0) is the specific heat 
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FIG. 19: Normalized deviation function calculated within the 
Eliashberg theory for an idealized two-peak phonon spectrum 
located at uji and 012 with equal electron-phonon coupling 
parameters Ai = A2 = 1 and strong total coupling parameter 
of Aph, tot = 2. Shown are results for uji/uj2 = 2,4 and 8. 
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FIG. 20: Normalized field-dependent Sommerfeld param- 
eter 7(_H')/7n plotted against iJ"///c2(0). Filled circles: 
7(iy)/7N = [cp(T,i/)-Cp(T,0)]/7N at T = 2K for dif- 
ferent applied magnetic fields. The black line is a fit of 

7(H)/7n = {H/H,2i0)f-^ using H,2{0) = 11 T and 7n at 
T = 2 K. 



in the Meissner state. Specific heat data for MgC^ gNia 
at T = 2 K were analyzed in order to derive the field 
dependence of ■y{H). In Fig. [201 the obtained j{H)/jj^ 
is plotted against H/H^^i^) using i?c2(0) = 11.0 T. 

The field data of Cp/T shown in Fi g. |2?11 can be de- 
scribed in accord with results from Ref. [l^ by the expres- 
sion 7/7N ~ {H/Hc2{0) f''^ which differs from the linear 
7(7?) law expected for isotropic s-wave superconductors 
in the dirty limit. 

A non-linear field dependence close to j{H) oc H^'^ 
has been reported for some unconventional superconduc- 
tors with gap nodes in the quasiparticle spectrum of the 
vortex state as YBa2Cu307r^ and in the heavy fermion 
superconductor UPta;^ but also in some clean s-wave 
superconductors as CeRu2r^ NbSe2^2ii£ and the boro- 
carbides RNisBsC (R = Y, Lu)!^!*^ Delocalized quasi- 
particle states around the vortex cores, similar as in d- 
wave superconductors, seem to be responsible for the 
non-linear j{H) dependence in the borocarbides42iffi 



D. The main superconducting and thermodynamic 
parameters 

In this subsection we collect the values of the main 
physical parameters we have found experimentally and 
compare them with available data in the literature. In 
order to make this comparison as complete as possible 
we estimate (calculate) , from our data and from those of 
Ref. 0, the lower critical field Hci(O) and the penetra- 
tion depth Al(0) at zero temperature adopting the appli- 
cability of the standard Ginzburg-Landau (GL) theory. 



TABLE II: Main superconducting and thermodynamic elec- 
tronic parameters for MgCNi3. 

Present work" Ref. Ref. 1^^ Ref. 



Tc [K] 


6.8 


6.4 


7.63 


7.3 


Hc2{0) [T] 


11 


11.5 


14.4 


16 


^^c(0) [T] 


0.18 


0.29±0.04 


0.19 


0.6 


Hci(O) [mT] 


11.3' 


23±7' 


10.0' 


12.6** 


eGL(O) [nm] 


5.47 


5.4 


4.6 


4.5 


k{0) 


43.3 


29.0±5.0 


54.0 


51 


Al(0) [nm] 


237.' 


154±26'' 


248.' 


230.0' 


Ac 

e& [K] 


31.4 
2.09 
292 


33.6 
1.97 
287 


30.1 
2.1 

284 




l-'ln [K] 


143 


135.= 


161.' 





"Using a parabolically extrapolated Hc2{0) value. 
'Using the artificial WHH estimate for _f/c2(0). 
'^Calculated. 



Measured. 
'Calculated (Eq. JHJ using Apj, : 



1.85, A,f = 0.43, At* = 0.13). 



Within this theory the penetration depth is given by the 
relation 



Al(0) =K(0)eGL(0), 



(17) 



where the Ginzburg-Landau coherence length ^gl (0) and 
the Ginzburg-Landau parameter k are related to the up- 
per and the thermodynamic critical fields as: 



&l(0) 
«(0) 



V$o/27ri/e2(0), 

gc2(0) 

V2i/c(0)' 



with the fiux quantum $o- With Hc2(0) = 11 T and 
iJc(O) = 180 mT (see Sec. ESI, Cgl(O) = 5.47 nm and 
k(0) = 43.3 are obtained. Using theses values in Eq. 
itTzjl . the penetration depth is estimated to be Al(0) = 
237 nm. Our calcula ted value agrees well with measure- 
ments performed bv IProzorov et alfe resulting A t,(0) = 
(250 ± 20) nm. It should be noted that lLin et alfe mea- 
sured a penetration depth of At, (0) = (128 . . . 180) nm 
for their sample (see also Ref. 1511) , possible consequences 
will be discussed in Sec. IVI Al To complete the criti- 
cal field analysis, the lower critical field i?ci(0) can be 
estimated using 

ffci(0)i?c2(0) = H^iO) (InK(O) + 0.08) . 

With Hc2{0) = 11 T and k(0) = 43.3 we get i?ci(0) = 
11.3 mT, agreeing, well with iJci(O) — 12.6 mT, mea- 
sured by I.Tin et aJl^ The results are shown in Tab. 
where for comparison result of Refs. and 152 have 

been included. Comparing these sets one finds a general 
qualitative accord. 
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VI. THEORETICAL ANALYSIS AND 
DISCUSSION 

Naturally, the obtained parameter set is model depen- 
dent. In this context even the case of relatively simple 
Fermi surfaces provides a difficult task to solve the full 
three(four)-dimensional EHashberg problem with given 

VF{k) and a'^F{ k , k' ,uj) for all physical quantities of in- 
terest. However, the solution of this problem can be suf- 
ficiently simplified for three practically important cases: 

1. the relatively simple standard isotropic single band 
model (ISB), where vp{k) is constant and the 
spectral function a^F depends only on the boson 
(phonon) frequency, 

2. a separable anisotropic single band model which 
exploits the so called first order Fermi surface har- 
monic approximation and 

3. the isotropic two-band model (ITB). The latter is 
a straightforward generalization of the ISB with re- 
spect to two order parameters. 

Due to the present lack of single crystal samples we will 
ignore the second issue. In addition, the cubic structure 
of MgCNi3 suggests only weak anisotropy effects. 



A. The isotropic single band analysis 

In the following section the electron-phonon coupling 
strength Aph is extracted from a simultaneous analysis of 
the upper critical field and the penetration depth in terms 
of the unknown impurity scattering rate 7imp [K] . Since 
the specific heat measurements do not clearly character- 
ize MgCNig as a one- or multi-band superconductor, the 
analysis starts within an ISB model. Within this model 
the upper critical field Hc2{0) is given hy^ 



H,2i0) [Tesla] = Hfr^{0) 



0.137i„,p [K] 



rc(l + Aph) 



where 



Hg{0) [Tesla] = 0.0237^i±|Bj^^^|pS 



(18) 



(19) 



and 7imp = vp/hmp is the scattering rate which de- 
termines the intrinsic resistivity {Ump denotes the cor- 
responding mean free path). The London penetration 
depth including the unknown impurity scattering rate 
7imp is given by an approximative formula 



Al(0) 



AL(0)^(l + ApO(^l + 0.7^j (20) 



valid for Aph < 2.5 (see App. El for the exact numerical 
expression) , with the bare clean Hmit London penetration 
depth 



Al(0) 



c 



197.3 



^pi [eV] 



(21) 



Using Eqs. lQH|l and (|2Dl, 7^ = 7o(l + Aph) and the 
experimentally determined quantities from Tab. lul we 
now will check the applicability of the ISB model. For 
this aim we consider the ratio 



6.77 X 10- 



■7n 



LmolK^ 



i/c2(0) [Tesla] • V 



A' 



(22) 

which includes the values of six experimentally read- 
ily accessible quantities: the Sommerfeld coefficient j^, 
Hc2{0), Tc, Al(0) and the volume of the unit cell. The 
dependence of R on the parameter 7imp/7c can be ex- 
pressed as 



R = 



1 + 0.357i„,p/A(0) 



(1 + Aph)" ' {1 + 0.137imp/ [Tc (1 + Aph)]} 



(23) 



In Fig. 121 the theoretical -R(7imp/7c) curves obtained 
from Eq. l|23ll for several Aph values are compared with 
the value of R derived from our experimental data which 
is represented in Fig. |21 as horizontal line. Crossing 
points between the theoretical prediction and the experi- 
mental result, which confirm the applicability of the ISB, 
are found for Aph > 0.8 at low scattering rates. Even in 
the case of higher electron-phonon coupling constants of 
Aph ~ 2, a clean limit scenario with 7imp/7c ~ 1 is fa- 
vored within the ISB analysis. Dirty limit (with weak or 
medium coupling) as proposed in Ref. |0can be excluded 
from the i?-check in Fig. [21 

It is noteworthy that the proposed i?-check is much 
more convenient than the similar Q-check, proposed re- 
cently by two of the present authors since the depen- 
dence on Aph is considerably weaker for R and, which is 
more important, R does not depend on the band struc- 
ture calculation. Thus comparing the results derived 
above with the expectations from these calculations, ad- 
ditional information on the nature of superconductivity 
in MgCNig, may be extracted. From Eqs. ifTsjl and ltT9)l . 
the effective Fermi velocity (in 10^ m/s) 



VF = 0.154 (l + Aph)^-^rc[K] X 




0.137imp [K] / [Tc [K] (1 -f Aph)] 



H,2{0) [Tesla] 



is obtained. Using the very weak scattering rates 
7imp/7c < 1 derived above and the experimental values 
Hc2{0) = 11 T and = 6.8 K, one estimates from Eq. 
l(24ll vf ~ (0.60 . . . 1.08) X 10 cm/s for electron-phonon 
coupHng constants in the range of 0.8 < Aph < 2.0. Com- 
paring this result with our band structure calculations 
(see Sec. InJ, one realizes consistence with the averaged 
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FIG. 21: Parameter R vs. 7imp/rc according to Eq. Ij23ll in 
the range of electron-phonon coupling constants 0.8 < Aph < 
2.0. Horizontal dash-dotted line: Experimental result for 
MgCi gNis derived from Eq. ll22ll . 
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FIG. 22: Left panel: Correlation between impurity scattering 
rate 7imp and electron-phonon coupling constant Aph.h derived 
from Eq. lO and Eq. ^ using -ffc2(0) = 11 T, Tc = 6.8 K, 
VF = 1.07 X lO'^ m/s, Al(0) = 237 nm, 2A(0) « 2.2 meV. The 
point of intersection of both curves marked by a filled square 
points to an electron-phonon coupling constant of Aph.h = 
1.74... 1.78 in the investigated MgCj gNis sample. Right 
panel: Ratio [[Hc2{0) / H^2{0)) — l] plotted against Aph,h. 
The filled square again corresponds to Aph,h = 1.74... 1.78. 
From ii'c2(0) = 11 T one estimates H^iO) « (8.79 . . . 9.07) T 
for the upper critical field in the clean limit. 



'i'hc2,h = 1-07 X 10'' m/s from the two hole Fermi surface 
sheets (Sec. InJ for strong electron-phonon coupling of 
Aph ~ 2.0. Thus, the relatively high value of the upper 
critical field of i?c2(0) = 11 T can be attributed to strong 
electron-phonon coupling for the hole subsystem. The 
second electron band plays a minor role for i?c2(0) due 
to its much faster Fermi velocities and the much lower 
partial density of states. 

Having adopted the dominant hole picture, we also can 
start from the band structure results, using the Fermi 
velocity wp.h and the plasma frequency Wpi.h of the hole 
band. Then we have to find consistent values of Aph and 
7imp, which describe the Hc2{0) and Al(0) data. 

From the plasma frequency of band 1, ftwpi i = 

1.89... 1.94 eV (see Sec. HD, we get Al(0) = 
(101.7. . .104.4) nm, using Eq. (EU. With the empir- 
ical values of Al(0) = 237 nm and 2A(0) « 2Aexp = 
2.2 meV = 25.5 K for the superconducting gap (see Sec. 
IV C|l . Eq. ll2(11l depends only on 7i„ip and Aph,h (of the 
hole band). The same applies to Eq. ifTsji . using the 
experimental values i?c2(0) — 11 T, = 6.8 K and 
the calculated average Fermi velocity of the hole band, 
^'hc2,h = 1.07 X lO*" m/s. The correlation between 7iinp 
and Aph,h, emerging from these two equations, is shown in 
the left panel of Fig. [221 The intersection of both graphs 
gives Aph,h = 1.74... 1.78 and 7inip = (31.0 ... 36.0) K. 
Thus, we arrive at a higher, more reaHstic scattering rate 
compared with 7imp ~ Tc obtained from the i?-check in 
Fig. ED The corresponding ratio [(i/c2(0)/i/c^^(0)) - l] , 
giving the deviation of i?c2(0) from the clean limit value 
H^l{0), is plotted in the right panel of Fig. [221 One gets 
H^llo) « (8.79... 9. 07) T. 

To summarize this part, already in the simplest possi- 



TABLE III: Characteristic phonon frequency and coupling 
parameters derived by analyzing the experimental data of the 
present MgC^ gNis sample. 







Cp analysis 
normal state si state 
(Sec. WM (Sec. [VC% 


Hc2 analysis 
(Sec. IVIAt 


^In 


[K] 


143 88 . . . 134 




Aph 




1.85 1.9... 2.3 


1.74. . . 1.78" 


Asf 




0.43 





"Limited to band "1". 



ble approach two general properties of MgCNi3 are de- 
rived: 

1. strong electron-phonon coupling and 

2. intrinsic clean limit at least for the hole subsystem. 

Nevertheless, it should be noted that recent prelimina ry 
measurements of the penetration depth by iLin et ahU ^ 
resulting in Al(0) = (128 . . . 180) nm are not compatible 
with the presented effective single band analysis (see also 
Ref. 113) ■ Especially the i?-check (Eqs. l|221l and ll23|l ) 
results in unphysical solutions (Aph = 30 as a lower limit), 
using the values presented in Ref. il^l and^l (see as well 
Tab. The consequences, if these measurements could 
be verified, remain unclear. 

B. Strong coupling and enhanced depairing 

Several results of our analysis of the experimental data 
are summarized in Table lllll The comparison of the esti- 
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mated Aph values clearly points to strong electron-phonon 
coupling. Nevertheless, the strong coupling scenario re- 
alized in MgCNig has been questionedii*i^ The strong 
electron-phonon coupling found for MgCNi3 requires a 
sizable depairing contribution to explain the low Tc value, 
otherwise at least a twice as large would be expected. 
It is illustrative to compare different approaches for the 
calculation of to analyze the electron-phonon coupling 
strength under consideration of the low temperature up- 
turn of the specific heat in the normal state (see Sec. 

In a first approach usually the low temperature De- 
bye approximation is used to extract the Debye temper- 
ature which we did in Sec. IIV Bl for comparison. Our 
result of = 292 K is in agreement with previous mea- 
su rements of iLin et alJ i^ lMa,o et alfe and calculations 
of llgnatov et al^ fit should be noted, that our specific 
heat data were corrected by carbon contribution (see Fig. 
EJ, without this correction we arrive at 8^ ~ 285 K). In 
this analysis the McMillan formula 



30 
25 
7 20 
3^ 15 
10 



tOp = 292 K 
(0|„ = 143 K 

tO|„ = 130 K , = 0.43 

tO|„ = 130 K , X„ = , n" = 0.4 




1.45 



exp 



-1.04 



1 + A 



ph 



Aph - (1 + 0.62Aph) 



(24) 



is usually appHed. This approximation is only useful for 
a special phonon spectrum with o^in/wD ~ 0.6. In the 
case of MgCNig we found « 0.30 . . .0.49 (corresponding 
to Win ~ 88 . . . 143 K) and the Allen-Dynes formula (Eq. 
(jHJ) should be applied instead. 

Fig. UHl compares both equations using 0^ = = 
292 K (dotted line) respectively win = 143 K (dash- 
dotted line). In both cases the Coulomb pseudopotential 
was fixed to ^* = 0.13. Apart from the deviation between 
Eq. © and Eq. (EU due to the ratio c^in/wo < 0.49, 
both equations seem to result in a moderate electron- 
phonon coupHng of Aph = 0.67. . .0.82 if no additional 
pair breaking effects are considered. 

However, we remind the reader, that the experimen- 
tal and theoretical picture of MgCNig strongly indicates 
strong electron-phonon coupling and a spin fluctuation 
contribution. The solid line compared to the dash-dotted 
line in Fig. |22l shows that the dependence of on 
Aph is strongly influenced by pair-breaking contributions 
such as the presence of enhanced electron-paramagnon 
coupling Agf = 0.43. The same situation in the imag- 
inable case of purely static pair-breaking, expressed by 
fx* = 0.4 is given by the dotted lin e. A very similar 
result was reported by llgnatov et alB ^ who proposed a 
phonon-softening scenario with Tc = 8 K, win = 120 K, 
Aph = 1-51 and an enhanced = 0.33 due to spin fluctu- 
ations. In any case the superconducting transition tem- 
perature is strongly suppressed by pair-breaking contri- 
butions. 

Tc of MgCNi3 would rise up to ~ 20 K, if one somehow 
could suppress the electron-paramagnon interaction. In 
that case the electron-phonon coupHng would not be af- 
fected and the dash-dotted line in Fig. l23lwould become 
reality. 



FIG. 23: Variation of Tc with Aph without electron- 
paramagnon interaction and "normal" Coulomb pseudopoten- 
tial fj,* = 0.13 according to Eq. 12411 (dotted line) and Eq. lO 
(dash-dotted line) and with enhanced pair-breaking contri- 
bution according to Eq. ^ hy n* = 0.13 and Agf — 0.43 
(solid line), respectively fi* = 0.4 and Agf = (dashed line). 
The characteristic phonon frequencies are chosen from Sec. 
IIV Bl fdotted line). Sec. IV Bl fdash-dotted line), respectively 
Sec. IV C\ (solid and dashed line). The points of intersection 
of the curves with the horizontal line at Tc = 6.8 K show 
the electron-phonon coupling strengths Aph resulting in the 
different approaches. 



Within the phonon-softening scenariojS^ which was re- 
cently observed in neutron-scattering measurements (21 a 
part of the low temperature speciflc heat anomaly may be 
of phonon origin (as stated in Sec. IVBII . In this picture 
the electron-paramagnon coupling would be reduced with 
the possibility of a paramagnon shift to higher tempera- 
ture s. This is consistent with uj\„ « 100 K (lower limit 
of the result from Sec. IV Cp and an electron-paramagnon 
coupling constant of Agf ~ 0.25. Using these numbers in 
Eq. (O the electron-phonon coupling constant amounts 
Aph ~ 1.6. To flnd the composition of the phonon and 
paramagnon contribution to the upturn, low temperature 
neutron-scattering measurements should be performed. 
In this context we remind the reader, that spin fluctu- 
ations are known to show a dependence on the applied 
magnetic fleld, which indeed is seen in Fig. 1121 



C. Multi-band effects beyond the standard 
approach 

Multi- (two-) band (and similar anisotropy) effects for 
several physical properties in the superconducting state 
are in principle well known for a long timepA especially 
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for weakly coupled superconductors in the clean limit. To 
the best of our knowledge their interplay with disorder 
and strong coupHng effects is less systematically studied. 
In particular this is caused by the increased number of 
input parameters and the necessity of a large amount of 
numerical calculations. 

The multi-band character in MgCNig is not only sup- 
ported by the band structure calculations, but also by 
experimental findings, and for instance, refiected by the 
large gap found in tunneling measurementsi^^ii and a 
smaller one seen in NMR measurementsii 

Like in MgB2 the effect of interband scattering is ex- 
pected to be weak due to the presence of well disjoint 
FSS. However, due to different contributions of the par- 
tial density of states compared to the case of MgBj, the 
two-band character of MgCNi3 is less pronounced. In 
Sec. IV CI the total electron-phonon coupling constant av- 
eraged over all Fermi surface sheets, Aph,tot was estimated 
by Eq. lfT?Hl at Aph,tot ~ 1-9- Considering the band struc- 
ture calculation presented in Sec. lUl this value is to be 
distributed among the two contributing bands according 
to 



X - X ^h(O) ^ , Afel(O) 



7V(0) 



N{0) 



(25) 



With Aph,h = 1-74 . . . 1.78 (see Sec. NTK^ . the coupling 
in the second band amounts Aph.ei — 2. 58... 2. 81. Ob- 
viously this strong mass enhancement in both bands is 
not compatible with the low value of Tc = 6.8 K. So, 
as in the single band case a pair breaking contribution is 
needed. 

Measurable quantities describing the superconducting 
transition in the case of MgCNi3 (particularly Eqs. Ijl4|l l 
are affected in opposite ways by strong coupling effects 
from one side and two-band effects from the other sides. 
Here we will briefly show, how these different effects in- 
fluence the speciflc heat jump Ac/7nT'c according to Eq. 
Ijl4b|l . Considering these effects. 



Ac 

is derived, where 



1.43F(^*) 



BiO^ (1 + W^^) 



(1 + v) (1 + riv/z^) 

(1 + Aph,el) Nel 

(1 + Aph,h) iVh 



(26) 



(27) 



contains electron-phonon coupling and multi-band cor- 
rections and 



Ael 

denotes the gap-ratio. The function F(/i*), given by 
1.15 [1 + 0.156 tanh(5Ai*)] 



(28) 



F(M*) 



1 4-0.156tanh(0.5) 



has been obt ained by analyzing numerical data derived 
by ICa,rbott(fe i under consideration of enhanced pair- 
breaking in terms of the Coulomb pseudopotential fi* 




FIG. 24: Left panel: Dependence of the specific heat jump 
Ac/7Nrc on the gap ratio Ah/Aei and the characteristic 
phonon frequency coin within the two-band description for 
fi* — 0.13 (Eq. I|26|l '). Right panel: Possible solutions for 
Ah/ Ael and ujin to reach the experimental specific heat jump 
for the two cases /i* = 0.13 and /j* = 0.4 (enhanced pair- 
breaking). The two vertical lines mark the range of expected 
Ah/ Ael values (see Sec. IVC|l . 



(valid up to /i* ~ 0.4). The different electron-phonon 
coupling constants in both bands may involve different 
characteristic phonon frequencies win. This situation is 
formally taken into account by the parameter rj, given by 



V 



Bi(i%^)/B/'^'"''^ 



The general result, depending on the gap ratio z and 
the characteristic phonon frequency win (for simpliflca- 
tion Win, el = ^\n.yi is assumed) is shown in the left 
panel of Fig. |21 (with Aph,h = 1-76, Aph^ei = 2.7 and 
F(m*) = !)■ The opposite effect of strong couphng from 
one side and two-band behavior from the other side is 
clearly seen. The right panel shows possible solutions 
for z and win to reach the experimental determined jump 
Ac/j^Tc = 2.09. Considering the results from the nor- 
mal state speciflc heat analysis (win ~ 143 K) and the su- 
perconducting speciflc heat analysis (1 < z < 0.8), best 
consistency is obtained with ^* w 0.4 in full agreement 
with the pair-breaking scenario. 

Now, the influence of two-band corrections on the pen- 
etration depth Al(0) should be checked, since the experi- 
mental value of Al(0) = 237 nm was ascribed to the hole 
band in Sec. IVI Al We start with the inverse squared to- 
tal penetration depth of a two-band superconductor. It 
is given as a sum of the two corresponding contributions 
from each band: 

Arot(0) = /hAi;^(0) + /eiA,f (0), 
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FIG. 25: Penetration depth at zero temperature derived from 
Eq. Ij29|l vs. scattering rate in the hole band using 7ei/7h ~ 
4.81 (see Eq. ((TJl). 



with the gap ratio dependent weighting factors 
/i(z, Ai, iV;), where i = h, el and /i = 1 in the case of 
equal gaps (i.e. z = 1). The corresponding specific 
plasma frequencies, c ouplin g constants, and gaps do enter 
each term (see Eqs. H2fll21|l l. Then the total penetration 
depth can be rewritten as 

AL,h(0) = Al(0) X 



1 + 



^el'^pl.el 



(1 + Ah)(l + ^^) 



Af,^^l.h(l + Ael)(l+^|^)J 



with 



1 + vjz^ 



(29) 



1 + vjz^' 



With the above determined values, 
(with V and z according to Eqs. 



we get L K. 0.88 

1). Fig. m 



shows the contribution of the hole band to Al (0) for the 
LDA calculation, using 7ei/7h ~ 4.81 (see Eq. Q). It 
is seen, that AL,h(0) exceeds the experimental value of 
-^l(O) = 237 nm by no more than « 18 % (in the case 
of 7h = 31.0 . . . 36.0 K - see Sec. I VI Ap . indicating only 
small influence of the electron band on Al(0). Neverthe- 
less its influence is not as small as in the case of -ffc2(0) 
and thus, we estimate an error of about ~ 10 % for the 
electron-phonon coupling constant of the hole band, de- 
termined in Sec. IVI Al 



VII. CONCLUSIONS 

Our analysis of MgCNi3 revealed a highly interesting 
interplay of different, on flrst glance unexpected adversed 



physical features or tendencies all present within one ma- 
terial causing a rather complex general behavior. This 
novel superconductor has been interpreted so far as stan- 
dard s-wave BCS superconductor or as unconventional 
superconductor with strong or medium electron-phonon 
coupHng. Last but not least, considerable pair-breaking 
contribution due to spin fluctuations and / or Coulomb 
repulsion have been suggested from theory and experi- 
ment. 

The present analysis is the flrst approach to reconcile 
the unusual experimental flndings within a unifled phys- 
ical picture. It reveals strong electron-phonon coupling 
combined with medium electron-paramagnon coupling. 
Strong electron-phonon coupling was derived from spe- 
ciflc heat data in the normal and superconducting state 
independently. An unusual upturn of the speciflc heat 
in the normal state observed at low temperatures can 
be attributed to spin fluctuations and / or a softening 
of low-frequency phonons. To specify the contribution of 
the belonging electron-boson interactions to the low tem- 
perature speciflc heat anomaly, low-temperature neutron 
measurements are necessary in order to investigate the 
evolution of the lattice excitations, which may even be 
modifled by the transition from the normal to the super- 
conducting state. 

The electronic speciflc heat data show an exponen- 
tial temperature dependence at low temperatures which 
is a strong indication for s-wave superconductivity in 
MgCNig. It was shown that a contribution of a second 
band could not be excluded but even complies with recent 
tunneling measurement results. The multi-band charac- 
ter of MgCNig is proved by band structure calculations. 
However, with respect to superconductivity the two-band 
character of MgCNig is much less pronounced than in 
the model compound MgB2. That is due to the predom- 
inance of a hole band with a large density of states in 
MgCNig, whereas in MgB2 the densities of states of both 
bands are comparable. Therefore, several properties such 
as the speciflc heat or the upper critical fleld can be de- 
scribed to flrst approximation reasonably well within an 
effective single band model. Nevertheless, other proper- 
ties such as the Hall conductivity and the thermopower 
require a multi-band description, i.e. at least one effective 
electron and one effective hole band (see App. El- Pre- 
vious theoretical analyses based on single-band models 
could describe only few physical properties. As a conse- 
quence of the oversimplifled approaches they blamed the 
local density approximation to fail seriously. This is in 
sharp contrast to our analysis of the upper critical fleld 
yielding an effective Fermi velocity agreeing well with the 
LDA hole band prediction. Our proposed effective strong 
coupling two-band approach explains the complex behav- 
ior observed for MgCNig and is expected to hold for other 
still not examined physical properties. 

The highly interesting interplay of strong electron- 
phonon coupling on multiple Fermi surface sheets, soft- 
ening of lattice excitations, the strong energy depen- 
dence of the density of states near the Fermi energy 
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of one band (van Hove singularity), and paramagnons 
or strong Coulomb repulsion for a realistic, anisotropic 
multi-band electronic structure with nesting features in 
this compound highly motivates further experimental 
studies. Investigating the influence of impurities or slight 
stoichiometry-deviations on the electronic and bosonic 
properties would be as helpful as making of purer sam- 
ples and single crystals to perform quantum oscillation 
studies like de Haas van Alphen measurements. 

Deepened theoretical studies are needed to clarify re- 
maining quantitative details and to extend the present- 
day strong coupling Eliashberg theory with all peculiar- 
ities of MgCNig. 



APPENDIX A: PENETRATION DEPTH : 
STRONG COUPLING AND IMPURITY 
SCATTERING 

We present a simple semi-analytic expression for the 
penetration depth at T = OK for type-H superconductors 
valid in the London limit. Thereby strong coupling and 
impurity scattering effects are treated on equal footing 
within standard Eliashberg theory. In calcul ating At, (0) 
we start from an expression proposed flrst bv lNamfe ^ ; 
later on frequently used in the literatur oi'?^n'?Tii?ff 



1 and 



(0) 



ttTlo, 



pi, 



(iWn) 



Z(itJ„)K2 -f-A2 (iWn)] 



3/2 ' 



(Al) 



where iw„ = i7r(2n — 1)T, n ~ 0, ±1, ±2, . . . are the Mat- 
subara frequencies and Z{iujn) as well as A(iw„) denote 
the renormalization factor and the gaps, respectively. 
The result of our numerical calculations of Eq. (IA1|1 com- 
pared with the approximation given by Eq. Ij20ll is shown 
in Fig. ESl 

One realizes only small deviations not exceeding 8 to 
10 % which is sufficient for our qualitative estimate of 
large mean free paths Zimp compared with the coherence 
length eGL(O). 



APPENDIX B: TWO-BAND INFLUENCE ON 
OTHER PHYSICAL QUANTITIES 

With the help of the two-band scenario even difficulties 
found explaining Hall data can be overcomei^ii^ Within 
the two-band model (Sec. the Hall constant is defined 
as 



i?H = 



((Tel -t- (Jhf 



using i?H,ei = — ^H,h with rih = riei (due to the even 
number of electrons per unit cell) and the ratio of the 
hole and electron conductivities 
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FIG. 26: Results of strong coupling calculations for the pen- 
etration depth at zero temperature (Eq. JAlIl l for several 
electron-phonon coupling constants A vs. impurity scatter- 
ing rate 7i mp (i n units of the gap Agxp = 1.1 meV as derived 
from Sec. IV C|l in comparison with the approximate expres- 
sion provided by Eq. 12011 . 



where 7h and 7ei are the corresponding scattering rates, 
we get 



1 



1 



(Bl) 



with i?H,ei = —1/ (jT-eiec). The number of charges per 
unit cell from LDA-FPLO calculations (see Sec. |n)| 
amounts n = 0.285 (comparable to n = 0.303 of Ref. 
0). The resulting theoretical charge carrier density of 
riei w 5.2 X lO^^/cm'^ and the calculated conductivity ra- 
tio of a; 1.403 (see Sec. should now be compared 
with measurements in terms of the effective Hall constant 
given by Eq. (jBTJ. 

Since so far reported samples are widely spread in 
terms of their residual resistivities, grain boundary ef- 
fects, affecting the Hall conductivity should be taken into 
account. This can be done in a first approximation by 
analyzing the mean free path (Eq. Q), considering the 
hole subsystem 



Up,h = 4.9 X 10^ 



VF,h [lO'' cm/s] 



(t^pLh [eV]) Pq{1 



[/iilcm] 



Similar considerations should be applied to the analysis 
of the thermopower, where also an effective electronic 
behavior has been observediifi 
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